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Abstract 

Using the X-ray data from ASCA, spectral and spatial properties of the intra- 
cluster medium (ICM) of the cD cluster Abell 1795 are studied, up to a radial distance 
of ~ 12' (~ 1.3 h'^Q kpc). The ICM temperature and abundance are spatially rather 
constant, although the cool emission component is reconfirmed in the central region. 
The azimuthally-avcragcd radial X-ray surface brightness profiles arc very similar be- 
tween soft (0.7-3 kcV) and hard (3-10 keV) energy bands, and neither can be fitted 
with a single-/? model due to a strong data excess within ^ 5' of the cluster center. 
In contrast, double-/? models can successfully reproduce the overall brightness profiles 
both in the soft and hard energy bands, as well as that derived with the ROSAT PSPC. 
Properties of the central excess brightness are very similar over the 0.2-10 keV energy 
range spanned by ROSAT and ASCA. Thus, the excess X-ray emission from the core 
region of this cluster is confirmed for the first time in hard X-rays above 3 keV. This 
indicates that the shape of the gravitational potential becomes deeper than the King- 
type one towards the cluster center. Radial profiles of the total gravitating matter, 
calculated using the double-/3 model, reveal an excess mass of ~ 3 x 10^^ Mq within 
~ 150 h^Q kpc of the cluster center. This suggests a hierarchy in the gravitational 
potential corresponding to the cD galaxy and the entire cluster. 



subject headings: galaxies: clusters: individual (Abell 1795) - galaxies: elliptical and lentic- 
ular, cD - X-rays: galaxies 



1 INTRODUCTION 



Imaging soft X-ray observations have established that nearly two thirds of clusters of galaxies, 
particularly those with cD galaxies, show "central excess X-ray emission" (Jones & Forman 
1984; Edge & Stewart 1991ab; Edge, Stewart, & Fabian 1992). That is, diffuse X-ray 
emission from their intracluster medium (ICM) exhibits centrally peaked surface brightness 
profiles when compared with /3 models (e.g., Sarazin 1988) which describe an isothermal 
gas sphere hydrostatically confined in a King-type potential (King 1972). Since the X-ray 
volume cmissivity from the ICM is proportional to the ICM density squared, the excess 
X-ray emission generally implies an excess ICM density in the cluster core region. This may 
in turn result from two apparently independent (and mutually non-exclusive) mechanisms. 
One is decrease in the average ICM temperature (or presence of a cooler ICM phase) in the 
cluster core region, particularly due to cooling fiows (see Fabian 1994 for a review); then 
the /3-model is no longer valid since the gas is not isothermal, and the central ICM density 
increases so as to maintain the pressure. The other is deviation of the gravitational potential 
from King-type ones, in such a way that the potential exhibits an additional drop in the 
cluster center which attracts an excess amount of ICM. 

Among these two mechanisms, the former effect is usually thought to be significant, 
because the ICM is often multi-phased in the central regions of cD clusters (Fabian et al. 
1994; Fukazawa et al. 1994; Matsumoto et al. 1996), and the radiative coohng time of the 
ICM in such regions is estimated to be generally shorter than the Hubble time (Fabian 1994). 
In contrast, the effect of potential shapes on the central X-ray excess has so far remained 
inconclusive, because the imaging X-ray telescopes were available only in energy ranges below 
~ 3 keV where the X-ray brightness is usually more sensitive to the temperature decreases 
than to the detailed potential shapes. 

With ASCA (Tanaka, Inoue & Holt 1994), the energy range available for X-ray imaging 
studies has been expanded up to 10 keV, and the spectroscopic capabilities have been drasti- 
cally improved. Observations of the Fornax cluster with ASCA have revealed that the strong 
central excess brightness around its cD galaxy, NGC 1399, is a consequence of the shape of 
the gravitational potential, which requires two distinct spatial scales corresponding to the 
entire cluster and the central galaxy (Makishima 1995, 1996; Ikebe 1995; Ikebe et al. 1996). 
The optical surface photometry of NGC 1399 shows essentially the same effect (Schombert 



1986). Similarly, the prominent central excess emission from the Hydra-A cluster, previously 
known in soft X-rays (David et al.l990), has been detected in harder X-rays up to 10 keV 
(Ikebe et al. 1997). Therefore the potential of this cluster is inferred to exhibit a marked 
central dimple to confine the excess hot ICM. Even in cD clusters with prominent central 
cool emission, such as the Centaurus cluster (Pukazawa et al. 1994; Ikebe 1995) and the 
Virgo cluster (Matsumoto et al. 1996), non-cooled ICM was found to fill a major fraction of 
the cluster core volume; two-phase modeling of the ICM in the Centaurus cluster has yielded 
a potential profile with a noticeable feature in the central ~ 50 kpc (Ikebe 1995). 

Similar effects have been observed from smaller systems as well. The giant eUiptical 
galaxy NGC 4636 was revealed to sit at the center of two concentric X-ray brightness compo- 
nents of similar temperatures, having distinct scales of ~ 30 kpc and > 200 kpc (Matsushita 
1997; Matsushita et al. 1996, 1997). Therefore NGC 4636 must be surrounded by hierarchi- 
cally nested two gravitating halos, one corresponding to the galaxy itself and the other to 
a galaxy group which is otherwise difficult to recognize. The nested two-component X-ray 
brightness distributions have also been detected with ROSAT from several galaxy groups 
(Mulchaey & Zabludoff 1997). 

Thus, recent X-ray observations suggest that the gravitational potential shape has a 
significant effect on the phenomenon of central excess X-ray emission. It may be that the 
cluster potential generally consists of doubly-nested hierarchical components, corresponding 
to the cluster and the cD galaxy. Alternatively, the potential shape of a self-gravitating 
system may normally exhibit a significant central cusp, as suggested by some gravitational 
lensing measurements (e.g., Wu & Hammer 1993; Miralda-Escude & Babul 1995) and N- 
body simulations (Navarro, Prenk, & White 1996). If either (or both) of these possibilities 
are true, the central excess brightness of clusters should be detectable in hard X-rays as a rule 
rather than the exception. In this paper, we examine this possibility using the prototypical 
central-excess cluster, Abell 1795 (hereafter A1795). 

A1795, a richness class 2 cluster at a redshift ol z — 0.0616, is one of the most equihbrated 
clusters with nearly circular X-ray isophotes (Briel & Henry 1996). The soft X-ray brightness 
within 4' — 5' of its cD galaxy (MCG -1-05-33-005) exhibits a strong excess, above a /9-model 
fitted to the outer-region brightness profile (Jones & Forman 1984; Briel & Henry 1996). 
In the 0.5-3 keV band, the central excess emission is so strong that it accounts for 30-50% 
of the total luminosity in the same energy band. Mainly based on these morphological 



results in soft X-rays, the cooling flow rate in this cluster was estimated to be as large as 
M ^ 478 h^^ Mq yr-i with EXOSAT (Edge, Stewart, & Fabian 1992), or M = 512 h^^ Mq 
yr~^ with ROSAT (Briel & Henry). Spectroscopically, the ICM of A1795 has an emission- 
weighted mean temperature of 5.3 keV (Hatsukade 1990; Tsuru 1992), together with a central 
cool emission component (Fabian et al. 1994; Mushotzky et al. 1995; Briel & Henry 1996). 
The spectroscopy with ASCA (Fabian et al. 1994) provided a somewhat lower cooling-flow 
rate of ~ 131 K^q Mq yr~^. Outer regions of A1795 are reported to be rather isothermal 
(Mushotzky 1994; Ohashi et al. 1997). Besides all these X-ray studies, an optical surface 
photometry of A1795 (Johnstone, Naylor, & Fabian 1991) also illustrates an elongated large- 
scale halo around the cD galaxy, which is similar to that of the Fornax cluster (Schombert 
1986). Clearly, this object is ideal to our search for the hard X-ray excess in the central 
region. 

Using the ASCA data of A1795, we have actually detected the central excess brightness 
in the hard energy band up to 10 keV, and found that the excess percentage in the 3-10 keV 
band is very similar to those found below ~ 3 keV. We also analyzed the archival ROSAT 
data and obtained a consistent result. Unless stated otherwise, we use the 90% confldcnce 
level throughout this paper, and assume the Hubble constant to be Hq — 50 h^Q km s~^ 
Mpc~^. We also employ the solar abundance ratios from Anders & Grevesse (1989), with 
Fe/H= 4.68 x 10"^ by number. 

2 OBSERVATION 

ASCA observations of A1795 were conducted for ~ 40 ks on 1993 June 16-17, during the 
PV (Performance Verification) phase. The field of view was centered at a = 13*^48^30^.4, 
5 = 26°37'39".7 (J2000). The data from the GIS (Gas Imaging Spectrometer; Ohashi et 
al. 1996; Makishima et al. 1996) were taken in the normal PH mode, and those of the SIS 
(Solid-state Imaging Spectrometer; Burke et al. 1994) were acquired in 4 CCD Bright/Faint 
mode. The same dataset was already utilized by Fabian et al. (1994), Mushotzky (1994), 
Mushotzky et al. (1995), and Ohashi et al. (1997). 

We screened the GIS and SIS events with the standard ASCA data selection criteria. 
Speciflcally, we required the telescope view direction above the Earth's limb to be > 10° for 
the GIS data, > 10° for the night-time SIS data, and > 20° for the day-time SIS data to 



avoid the light leakage. The minimum cutoff rigidity was set to be 8.0 GcV c^^ for both the 
GIS and the SIS. These screening processes have provided ~ 35.8 ks and ~ 37.3 ks of good 
exposures with the GIS and the SIS, respectively. 

Figure 1 shows GIS images of A1795 in the 0.7-3 keV and 3-10 keV bands, where 
the cosmic X-ray background and internal detector background have been subtracted by 
utilizing observations onto several blank sky fields. The two images arc very similar to each 
other, both having a good circular symmetry. The X-ray brightness centroid is found at 
a = 13^48"'48^9 and 5 = 26°36'20".4 (J2000), with a typical uncertainty of 1'. This agrees 
with the optical position of MCG -1-05-33-005 within ~ 1.1'. The X-ray emission extends 
at least to 20': considering that 1' corresponds to ~ 108 h'^Q kpc, the ICM emission is thus 
detectable at least up to 2.2 h^Q Mpc. 

A Seyfert 1 galaxy, EXO 1346.2-^2645, is located about 6' south of the A1795 center and 
is included in the field of view of both the GIS and SIS. Since the X-ray flux of this object 
is about 6% of that of A1795, we exclude a region around it in the following data analysis. 
There is yet another active galactic nucleus (AGN) to ~ 14' south of A1795, emitting a 
flux which is 2% of A1795. We accordingly limit our analysis to within 12' from the cluster 
centroid. Although the extended point spread function (PSF) of the ASCA X-ray telescope 
(XRT; Serlemitsos et al. 1995; Tsusaka et al. 1995) causes some fraction of X-rays from 
these AGNs to leak into the analysis region, the effect is confirmed to be neghgible compared 
with the photon counting errors. 

3 SPECTRAL ANALYSIS 
3.1 Extraction of the Spectra 

Using both the GIS and SIS, we accumulated photons in five annular regions of r < 2', 
2' < r < 4', 4' < r < 6', 6' < r < 8', and 8' < r < 12', where r is the projected radius 
from the X-ray brightness centroid. However, due to the smaller field of view of the SIS 
(~ 22' X 22') and an offset of the cluster centroid from the detector center (Fig.l), the SIS 
data only partially covers the regions outside r — & . The background was subtracted from 
each spectrum, in a similar way as mentioned in §2. Corresponding spectra from the two 
GIS detectors (GIS2 and GIS3) were then combined together; so were those from the two 
SIS detectors (SISO and SISl). We show these spectra in Figure 2, where we combined four 



outer-region spectra for presentation. 

As seen in Figure 2, all these spectra exhibit strong Fe-K emission lines. In addition, 
both instruments reveal somewhat enhanced emission at ~ 1 keV in the spectra from the 
inner region (r < 2') , which is indicative of an additional cool component as already reported 
by Fabian et al. (1994). However in the outer regions (r > 4'), the spectral shapes do not 
depend significantly on the radius. Thus, the ICM of this cluster has relatively uniform 
properties in outer regions, in agreement with reports by Mushotzky (1994) and Ohashi et 
al. (1997), who studied temperature and abundance distributions using the same ASCA 
data. 

3.2 Conventional Spectral Fitting 

It is well known that the extended PSF of the ASCA XRT causes difficulties in the spectral 
analysis of extended sources (e.g. Takahashi et al. 1995; Markevitch et al. 1996; Ikebe 
et al. 1997; Ohashi et al. 1997). In the case of clusters of galaxies, X-rays from bright 
central regions are scattered by the PSF into fainter peripheral regions of the target, thus 
contaminating the outer-region spectra. Furthermore, the scattered photons form a harder 
spectrum than their original one, because the PSF has more prominent outskirts in higher 
energies. This often produces an artificial outward temperature increase. 

In spite of these problems, we first performed a brief spectral analysis in a conventional 
way just to grasp gross properties of the X-ray emission. That is, we jointly fitted the SIS 
and GIS spectra at each radius with the single temperature Raymond-Smith emission model 
(Raymond & Smith 1977), while the PSF effects are not taken into account. The absorption 
column density and the heavy-element abundances were left free, whereas the abundance 
ratios were constrained to obey the solar ratios. 

This conventional analysis yielded mostly acceptable fits to the spectra, and the tem- 
perature turned out to be in the range 5-8 keV with a mild outward increase due to the 
instrumental artifact (Takahashi et al. 1995). Also there is a shght temperature drop at the 
center, due probably to the cool component (Fabian et al. 1994; Briel & Henry 1996). Ex- 
cept these, the cluster appears to be rather isothermal in agreement with Mushotzky (1994) 
and Ohashi et al. (1997) who used the same ASCA dataset. The metaUicity is consistent 
with being radially constant at ~ 0.3 solar. There is some hint of excess absorption in the 
inner regions up to a few times 10^° cm~^, above the Galactic line-of-sight value of 1.1 x 10^° 



cm ^, but the effect is not as significant as in some other clusters (White et al. 1991). 

^From these fits, we estimate the X-ray luminosity of A1795 to be (10.0 ±0.1) x 10^^ h^Q 
erg s-^ in 0.5-3 keV or (10.4 ± 0.1) x lO^'' h^^ erg s"^ in 2-10 keV, both integrated up to 
r = 12' (1.3 h^Q Mpc) and corrected for absorption. The former value is close to the 0.5-3 
keV luminosity of (7.25 ± 0.07) x 10^^ erg s~^ obtained with Einstein within r = 12' (Jones 
& Forman 1984). Similarly, our 2-10 keV luminosity is in a good agreement with the value 
of 1.1 X 10^^ erg s~^ (Hatsukade 1989), obtained with the Ginga LAC in the same 2-10 
keV range but over a much wider field of view of 1° x 2°. Therefore, contribution to the 
X-ray luminosity from regions outside r ~ 12' is inferred to be negligible. Note that these 
Einstein and Ginga measurements both refer to h^o — 1, and that these ASCA luminosities 
are approximately free from the PSF effects. 

3.3 Comparison with a Simulated Isothermal Cluster 

In order to derive the ICM temperature distributions considering the PSF effects, we must 
solve the spectral couplings among different regions on the focal plane (e.g. Markevitch et al. 
1996; Ikebe et al. 1997). However such a detailed spectral investigation is beyond the scope 
of this paper. Instead, we examine the ASCA spectra against the hypothesis of constant 
ICM temperature, taking the PSF effects fully into account. 

For this purpose, we performed a Monte-Carlo simulation (Takahashi et al. 1995; Ikebe 
1995) which converts a model cluster into the simulated ASCA data. That is, we generate a 
number of Monte-Carlo photons by specifying their spatial and spectral probability distribu- 
tions. The photons are then "reflected" by the XRT and "detected" either by the CIS or the 
SIS, according to the known instrumental responses. This tells us how a "known" cluster 
looks like when observed with ASCA. The XRT-I-GIS angular response is represented by 
actual images of Cygnus X-1 (a bright Galactic X-ray binary) obtained with the CIS (Taka- 
hashi et al. 1995). However Cygnus X-1 is too bright for the SIS, so that this Monte-Carlo 
method cannot be apphed at present to the SIS data. 

In the simulation, we employ the Raymond-Smith model to express the spectral proba- 
bility distribution of Monte-Carlo photons. We assume the model cluster to have a free but 
spatially constant temperature, and flx the model cluster abundance at 0.30 solar according 
to the results of §3.2. The spatial probabihty distribution of Monte-Carlo photons is modeled 



by a /3-model, which describes the radial X-ray surface brightness profile as 

(1) 

where S{0) is the central brightness, r is again the projected radius, a is the core radius, and 
/? is so-called f3 parameter. We chose a = 0'.8 and j3 = 0.56 for our simulation. Later we 
in fact find that the brightness profiles of A1795 cannot be described by a single-/3 model. 
However, the /3-model employed here can describe the ASCA and ROSAT brightness profiles 
of A1795 to within 5% as we show later. The difference between this approximation and the 
true brightness distribution produces negligible effects in the simulated spectra. 

We generated about 10^ simulated GIS events, which were subsequently sorted into 5 
annular-region spectra exactly in the same way as the actual data. Then, at each radius 
range, we directly compared the actual and simulated GIS spectra, by subtracting the latter 
(after an appropriate re-scaling) from the former in channel-by-channel way, and evaluating 
residuals via the chi-square tests. The re-scaling factor was adjusted at each radius to 
minimize the chi-square. We then varied the model cluster temperature, to optimize the 
overall fit goodness. The model cluster with a temperature of 6.2 keV has been found 
to best reproduced the GIS spectra, yielding xV^ = 82/106, 92/106, 65/106, 53/106 and 
82/106, from the inner to outer annular regions. Therefore we conclude that the observed 
GIS spectra are consistent with those expected from an isothermal cluster at a temperature 
of 6.2 keV, when the PSF effects are fully taken into account. This "best-fit" temperature 
is estimated to be uncertain by ~ 0.3 keV. 

Figure 3 shows the ratio (instead of difference) between the actual and simulated GIS 
spectra, in the innermost and 6' < r < 8' radius ranges. As expected from the above results, 
the ratios are approximately constant across the energy range. 

3.4 Central Cool Component 

Although the overall cluster is nearly isothermal as confirmed above, the central cool com- 
ponent detected by Fabian et al. (1994) reveals itself in the form of enhanced Fe-L complex 
(at ~ 1 keV) at about 1 keV, in the innermost spectra of Figure 2 and the spectral ratios 
of Figure 3a. To investigate this cool component, we extracted the SIS and GIS spectra 
from the inner region of r < 3'. This region, having by far the highest brightness, is little 
(< 5%) affected by photons which are scattered from outer regions by the PSF into this 
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region. (However, outer regions are strongly affected by photons from the central region.) 
We are therefore justified to analyze the spectra using ordinary energy responses calculated 
at the position of the cluster centroid, without considering X-rays from outer regions due to 
the PSF effects. 

We tried to fit these spectra jointly with a Raymond-Smith model whose temperature is 
fixed at 6.2 keV, according to the result of §3.3. However as shown in Table 1, the fit was 
totally unsuccessful even when the absorbing column density and the overall metallicity were 
allowed to float. We then let the temperature to vary freely. This gave a much improved 
fit together with a somewhat lower temperature (Table 1), which is close to that measured 
with Ginga (5.3 keV), but the fit still remained unacceptable. We therefore introduced 
another Raymond-Smith component of lower temperature, which is constrained to have the 
same abundance and the same absorption as the first component. This two-temperature 
(2T) fit has been mostly successful as shown in Table 1. The hot-component temperature 
(6.54lo;5o keV), thus derived jointly with the two instruments, agrees with that obtained 
via the simulation in §3.3. We also tried another two-temperature modeling after Fabian et 
al. (1994), who assumed that the hot component is absorbed with the Galactic line-of-sight 
column whereas the cool component is further absorbed by an excess column density. This 
model gave as good a fit {'X^ /v — 380/275) to the data, and the derived results are essentially 
the same as obtained above. 

/,From the two- temperature fit, we calculate the 0.5-3 keV luminosity of the cool com- 
ponent to be Lcoo\ = (1-43 ± 0.39) x lO'*'' h^Q erg s"'^, after removing absorption. This is 
14 ± 4% of the total 0.5-3 keV luminosity described in §3.2. The alternative modehng used 
by Fabian et al. (1994) yields a consistent value; Lcooi = (1.62 ± 0.46) x 10^^ h^Q erg s~^ 
also after removing absorption. When these values of I/cooi are substituted into equation 
(2) of the review by Fabian (1994), we get mass deposition rates of ~ 100 Mq yr~^, which 
are consistent with those obtained via the direct spectral fitting in terms of the cooling-flow 
model. We cannot however compare our cool-component luminosities with that derived by 
Fabian et al. (1994), which is not given explicitly in their publication. 

In order to more directly confirm the consistency with the result of Fabian et al. (1994), 
we tried an isothermal model plus an absorbed cooling flow model on the central (r < 2') 
SISO spectrum. The mass deposition rate turned out to be I38I52 h^Q Mq yr^\ with 
X^/i/ = 306/230. This reconfirms the report by Fabian et al. (1994) who used the same SIS 



data. 



4 RADIAL BRIGHTNESS PROFILES 

4.1 X-ray Brightness Profiles 

In Figure 4, we show radial count-rate profiles of A1795 in two representative energy bands, 
without correction for the XRT vignetting or the XRT PSF effects. They were derived by 
azimuthally averaging the counts from GIS2+GIS3 and SISO separately around the X-ray 
centroid, and then subtracting the corresponding backgrounds. In deriving these profiles, 
we excluded a region of radius 2'. 5 centered on the Seyfert 1 galaxy EXO 1346.2-1-2645. For 
comparison, we also show the instrumental PSF of the XRT-I-GIS combination. Although 
this PSF was calculated over the total 0.7-10 keV band, it can well approximate those in 
either soft or hard energy band. 

In order to compare the classical soft-band brightness profiles of A1795 against the hard- 
band ones that have been obtained for the first time, we have normalized the latter to 
the former at each radius. The results, or radial profiles of the spectral hardness ratio, 
are presented in Figure 4c for the GIS and SIS separately. They are very flat; the GIS 
profile is almost constant within ~ 10%, except beyond ~ 8' where the instrumental artifact 
mentioned in §3.2 begins to affect the data. The SIS hardness ratio indicates a slight drop 
in r < 2'; this is presumably because the SIS has a higher spectral sensitivity and a better 
spatial resolution to detect the central cool emission, whose presence was already confirmed 
in §3.4. The constancy of the hardness-ratio profiles implies that the central excess emission 
that has been observed in the soft X-ray band should exist also in the hard X-ray band. 
This inference is consistent with the results obtained in §3 that the spectral changes at the 
cluster center are rather modest. 

Two obvious criticisms may be addressed to the above inference. First, the central excess 
emission may have too soft a spectrum to be detected with ASCA. This is however unlikely, 
since the average ICM temperature of A1795 derived with ROSAT through the cluster center, 
2.9 keV (Briel & Henry 1991), is clearly high enough for the ASCA bandpass. Alternatively, 
the central soft X-ray excess may be undetectable with ASCA because of the heavy smearing 
by the PSF. In order to examine this issue in the simplest manner, an appropriately scaled 
3-10 keV radial profile of a point source was added to, or subtracted from, the hard-band 



radial profile of A1795. This is because any central excess emission should appear almost 
pointlike to ASCA. For the point source, we used the actual GIS data of Cygnus X-1 observed 
at the same position as A1795 (Takahashi et al. 1995). The hard-band profile of A1795, 
thus modified, was then divided again by the soft-band GIS profile of A1795. The results 
shown in Figure 4d clearly confirm that the hardness ratios deviate significantly from the 
actually measured values when the point-source contribution (either positive or negative) 
exceeds ~ 3% of the A1795 counts integrated within 12'. 

^From these considerations, we conclude that the hardness-ratio profile obtained with 
ASCA is very sensitive to small differences in the fractional central brightness excess between 
the two energy bands. Given that the relative excess is 30-50% in energies below ~ 3 keV 
(Jones & Forman 1984; Briel & Henry 1996), the relative excess should be at least 25% in 
hard X-rays. 

4.2 Single-/? Fitting 

4.2.1 The Entire Cluster Region 

To strengthen the conclusion derived in §4.1, we need to quantify the X-ray count-rate pro- 
files. However, an X-ray image obtained with ASCA is a very complex integral transform 
(not a simple convolution) of the original surface brightness distribution on the sky plane, 
because the PSF of XRT is axially non-symmetric and heavily position dependent. Ac- 
cordingly, instead of trying to inversely solve this integral transform, we again employ the 
Monte-Carlo simulation used in §3.3. That is, we generate simulated cluster data according 
to the j3 model of equation (1), in which the XRT-I-GIS PSFs are again represented by the 
actual GIS image of Cygnus X-1 observed at various detector positions (Takahashi et al. 
1995). We then derive simulated GIS count-rate profiles in several energy bands, in the 
same way as has been done with the actual GIS data. By evaluating their agreements with 
the actual GIS radial profiles, we search for the best-fit (3 model. We added a systematic 
error of 3% to the simulated radial profiles. This technique has been developed and applied 
to the ASCA data by Takahashi et al. (1995), Ikebe (1995), Markevitch et al. (1996), Honda 
et al. (1996), Matsuzawa et al. (1996), Ohashi et al. (1997), and Ikebe et al. (1997). 

We first apphed this /3-model to the GIS radial count-rate profiles over the whole radius 
range of r < 12'. As summarized in Table 2, rather flat /3 models with small core radii 



provided the best-fit solutions in the soft (0.7-3 keV) and hard (3-10 keV) energy bands. 
However, none of these fits are acceptable, exhibiting subtle but statistically significant fit 
residuals as shown in Figure 5a-5b. These failures of the fit can be ascribed obviously to the 
presence of the central excess emission, as already pointed out by Jones & Forman (1984) 
and Briel & Henry (1996). 

We also analyzed the processed archival ROSAT PSPC data of A1795, provided by the 
ROSAT Guest Observer Facility at the NASA/GSFC. These data, acquired on 1991 July 1 
and 2, comprise a part of the datasets used by Henry & Briel (1996). The pointing position is 
a = 13M8'"52.8^ 5 = +26°37'12".0 (J2000), just near the center of the cluster. By using the 
exposure map, we corrected the broad-band (0.2-2 keV) image for vignetting and variations 
in the exposure across the field of view. We then derived the azimuthally-averaged radial 
count-rate profile in the 0.2-2 keV range, in which the Seyfert 1 galaxy (EXO 1346.2-1-2645) 
was omitted by a circular mask of 1.'5 radius. The background is not subtracted, but 
considered in the latter model fitting. 

We performed a single-/? model fitting to the derived ROSAT count-rate profile, over 
the radius range of r < 12'. We approximated the PSF with a delta-function, which is 
justifiable when we ignore spatial structures much finer than ~ 1' (H. Bohringer, a private 
communication). We included a background component into the fit model assuming it to 
be free but spatially constant. As shown in Table 2 and Figure 5c, the fit turned out to be 
unacceptable, confirming the report by Briel & Henry (1996). 

In spite of these failures, we note that the single-/? model can approximate the observed 
X-ray profiles to a reasonable accuracy (e.g. within 5-10 %); this allowed us to utihze 
the single-/? approximation of the brightness profile in running the Monte-Carlo spectral 
simulation (§3.3). It is still more important to note that the best-fit /?-model parameters 
(though unacceptable) are quite similar between the ROSAT PSPC and the ASCA GIS, 
and among the two energy bands of ASCA (Table 2). Furthermore, the fit residuals behave 
essentially in the same manner in the three panels of Figure 5: the data exceed the model at 
r = 3' — 6', while the opposite occurs in other regions. These results support our inference 
made in the previous subsection that the count-rate profiles do not depend on the energy 
significantly. 



4.2.2 Outer Cluster Region 

When the central cluster region of ~ 5' radius is ignored, the single-/? model has given 
acceptable fits to the GIS profiles in the three energy bands (soft, hard, and total), as well 
as to the 0.2-2 keV ROSAT profile. We present these fits in Figure 6, and summarize the 
best-fit parameters in Table 3. 

Clearly, the /3-model parameters derived with the ASCA GIS are again nearly energy- 
independent within the respective errors. They also agree with those derived with the 
ROSAT PSPC which in turn are close to those obtained by Briel & Henry (1996) who used 
a larger ROSAT datascts, discarding the central ~ 5'. Furthermore, the previous results 
from the Einstein IPC (Jones & Forman 1984) are roughly consistent with those obtained 
here, as we can see in Table 3. In a word, regardless of the used instrument or the employed 
energy range, we obtain a core radius ~ 3' — 4' and a slope oi P — 0.8 — 0.9 for A1795, as 
long as the central ~ 5' region is ignored. 

When extrapolated inwards, these best-fit /? models fall significantly below the observed 
brightness (Fig. 6), thus reconfirming the prominent central excess emission. The last 
column of Table 3 gives a rough measure of this central excess, integrated within r < 5' 
and normalized to the total counts integrated up to a projected radius of 12'. Thus, the 
measurements with Einstein, ROSAT, and ASCA consistently indicate a central excess 
emission that reaches 20 — 30% of the total X-ray emission within 12', although the excess 
may decrease mildly with energy, which is most likely to be caused by the cool emission 
component. 

Combining these results with those obtained in the previous subsections, we draw the 
following two conclusions of basic importance. One is that the central excess emission of 
A1795 it too strong to be smeared out by the PSF of ASCA, and the other is that the central 
X-ray excess emission of this cluster is significantly seen over the entire 0.2-10 keV energy 
range spanned by ROSAT and ASCA. 

4.3 Double-/? Fitting 

In order to describe the central excess emission more quantitatively, we further attempted 
to fit the count-rate profiles with a sum of two P model components, a more extended one 
representing the emission outside ~ 5', and a compact one representing the central excess. 



We have found that this double-/? model gives acceptable fits to the profiles obtained with 
the ROSAT PSPC and the ASCA GIS, as summarized in Table 4 and shown in Figure 
7. All the observed count-rate profiles have thus yielded very similar parameter values for 
the two /3-model components. The more extended component is similar to that obtained in 
the single-/? fits excluding the central region (Table 3), while the compact component has a 
core radius which is close to that obtained with the single-/? fits including the central region 
(Table 2). 

We remark that the very central (r < 1') region of the cluster may remain unresolved even 
with the ROSAT PSPC. Consequently, the core radius of a ~ 1' for the compact component 
may be subject to instrumental angular resolutions, and the finer-resolution X-ray brightness 
profiles, e.g. those obtained with the ROSAT HRI, may no longer be adequately described 
by our double-/? model. Nevertheless, our principal conclusions should remain unchanged, 
because they all deal with angular scales of ~ 1' or larger. 

In Table 5, we summarize the luminosity, I/excess, carried by the narrower /? component. 
It amounts to 40 — 50% of the total (i.e., narrow plus wide) luminosity, Ltotab in both the soft 
and hard energy bands. This is refiected in the fact that the two /? components crossing over 
in Figure 7 have similar normalization ratios in the two ASCA energy bands. The excess 
luminosity defined in this way is systematically larger than those given in Table 3, by up to 
a factor of ~ 2. However, this is simply the matter of how to define the "excess". 

4.4 Fitting with the Universal Halo Model 

Through extensive A'"-body numerical experiments of structure formation in the cold dark 

matter (CDM) scenario, Navarro, Frenk, & White (1996) reported that the computed CDM 
mass density pdm takes a "universal halo shape" of 

where S and Eg are parameters. Makino, Sasaki, & Suto (1997) showed that the density of 
an isothermal gas hydrostatically confined in the potential of equation (2) can be expressed 
analytically as 

n{R) = no (l + 

where Uq and B are also parameters. By numerically projecting the square of this n{R) onto 
2 dimensions, we can obtain the X-ray surface brightness, 5'uh(^), emitted by an isothermal 



R\ R 



gas sphere confined in a universal CDM halo. 

Since the CDM distribution of equation (2) is more centrally peaked than the King- 
type distribution, the idea of universal CDM halo may explain the central excess of A1795. 
In order to examine this possibility, we fitted this -S'uh('^) to the ROSAT PSPC profile. 
However as shown in Figure 7d, the fit turned out to be as poor as that with the single- 
13 model (Fig. 5c), with y^/v = 56/16. The GIS profiles could not be fitted with this new 
model, either. This result is not surprising, since Syuir) is in fact fairly close to the ordinary 
/3-model profile S{r) given by equation (1), with only moderate excess towards the center 
(Makino, Sasaki, & Suto 1997). In short, the central excess emission from A1795 is too 
strong to be explained by the universal halo model as long as the isothermahty holds. 

5 DISCUSSION 

5.1 Summary of the Results 

Through a series of analysis of the radial count-rate profiles of A1795 obtained with the 
ASCA GIS and the ROSAT PSPC, we have confirmed that the central excess brightness of 
this cluster, which was well known previously in the soft X-ray band, does exist also in the 
3-10 keV band with a similar fraction. If employing /9-models, two distinct components are 
required to express the projected radial count-rate profiles obtained with the ASCA GIS in 
the energy range both below and above 3 keV, as well as that obtained with the ROSAT 
PSPC. 

We also carried out spectroscopic studies of the ICM in A1795 using the ASCA data. In 
essence, the ICM in A1795 is isothermal at a temperature of 6.2 kcV, except in the innermost 
region of radius ~ 3' (~ 320 h^Q kpc) where an additional cool emission component is seen 
as already report by Fabian et al. (1994). Employing the two-temperature fit obtained in 
§3.4, we calculated the luminosities of the cool component in various energy bands, as listed 
in Table 5. 

In the 0.5-3 keV band, these morphological and spectroscopic results compare as Lcooi/Lexc 
0.42 ± 0.16 (according to Table 3), or 0.24 ± 0.07 (according to Table 5). These ratios may 
not necessarily imply a discrepancy between L^ooi and I/excess > since our estimation of Lcooi 
depends on a particular modeling (i.e., two-temperature assumption) of the ICM. However 
the discrepancy becomes much more significant when the comparison is made in the newly 



explored 3-10 keV band; Lcooi/-^^excess = 0.15 ± 0.04 in reference to Table 3, or 0.07 ± 0.02 
in reference to Table 5. These small ratios directly result from the fact that the spectral 
cool component has a temperature of at most ~ 2 keV (Table 1), whereas the morphological 
excess component has a temperature rather close to the cluster-average value of 6.2 keV. 

Our results consistently indicate that the whole cluster volume of A 1795 is permeated 
with the hot isothermal ICM, including the central region where a small amount of cool 
component coexists. This is similar to the case of several other cD clusters (Fukazawa et al. 
1994; Ikebe 1995; Matsumoto et al. 1996; Ikebe et al. 1997). The central excess brightness of 
A1795 is inferred to be produced primarily by a central excess in the density, hence pressure, 
of this hot ICM itself, with the cool emission component giving subsidiary effects. To confine 
this excess ICM pressure, the gravitational potential in A1795 must exhibit an additional 
central drop as compared to the King-type potential. As indicated by Figure 7d, this central 
drop may be much deeper than is predicted by the universal halo model of equation (2), as 
long as the isothermality approximation is valid. 

5.2 Integrated Mass Profiles 

Referring to the double-/? solutions which provide acceptable fits to the observed brightness 
profiles, we quantify the shape of the gravitational potential. Specifically, we employ the 
best- fit double- /5 model obtained in the CIS 0.7-10 kcV band (Table 4), and back-project 
the constituent /9-model components individually into 3-dimensions. This provides emissivity 
profiles of the extended and compact components, denoted as ee{R) and e^R), respectively. 
These can be converted to the ICM density profile n{R) through the relation 



where it! is the 3 dimensional radius and A(T, A) is the cooling function for a given temper- 
ature T and abundance A. We then calculate the radially integrated total gravitating mass, 
Mtot(-R), using the usual equation of 



The total mass profiles calculated via equation (4) are presented in Figure 8, where we 
assumed two alternative cases for the temperature profiles T{R); either a constant ICM 
temperature at 6.2 keV, or a centrally decreasing profile of the form of 



ee{R)+e,{R) = A{T,A)n{Rf , 



(3) 





T{R) = 4.50 + OMR - 1.57 exp [-0.67i?] , 



(5) 



which approximately reproduces the ROSAT PSPC result of Briel and Henry (1996). Here 
the 3 dimensional radius R is in arcminute and the temperature T{R) is in keV. We also 
present the ICM mass profile, Micm(-R), calculated by integrating n{R) of equation (3). 

Inside ~ 50 — 100 h^^ kpc, the mass curves in Figure 8 should be regarded as subject 
to the limited angular resolution. Finer details of the central density profile could be studied 
using the ROSAT HRI data, but the lack of associated spectral information hampers the 
calculation of gravitational potential there. In this paper, we therefore limit our study to 
the information available with the ROSAT PSPC and ASCA. 

Figure 8 reveals a shoulder-hke structure in the Mtot(-R) profile at Rx ~ ISO/igo^ kpc, 
regardless of the assumed temperature profile. As a result, the total gravitating mass within 
Rx apparently exhibits an excess by Mcxcess — 3x 10^^ Mq, above the large-scale smooth mass 
distribution carried by the extended /9-component. This Mexcess is thought to be responsible 
for the central dimple in the gravitational potential. 

5.3 Nature of the Central Excess 

The shoulder-like structure seen in the Mtotai curve of Figure 8 may allow two shghtly 
different interpretations. One possibility is that the overall mass distribution is made up 
of a single mass component, presumably dominated by dark matter, whose density keeps 
rising to the center instead of forming a flat core. Although equation (2), representing 
one particular example of such cases, failed to explain the X-ray data under the isothermal 
assumption, some other numerical simulations (Moor et al. 1997) suggest even a steeper 
CDM halo profile. When taking into account the effect of the cool component and baryonic 
mass associated with the cD galaxy, such a centrally-steepening CDM halo may still be 
consistent with the observation. 

Alternatively, the Mtot(-R) profile in Figure 8 may be interpreted as exhibiting a gravi- 
tational hierarchy between the overall cluster and the cD galaxy. Such a nested potential 
structure was suggested or assumed previously (e.g. Thomas et al. 1987), and has been 
confirmed in the Fornax cluster (Makishima 1995; Ikcbc et al. 1996), the elliptical galaxy 
NGC 4636 (Matsushita 1997; Matsushita et al. 1996, 1997), and possibly in about 8 galaxy 
groups (Mulchaey & Zabludoff 1997). Then, the region inside Rx is considered to belong 
to the cD galaxy, while the outside region to the entire cluster. In reference to the above 
Mexcess and the optical B-band luminosity of the cD galaxy (6.5 x 10^^ in unit of the solar 



B-band luminosity), the mass-to-light ratio within Rx becomes ~ 50. 

The hierarchical potential structure is also suggested by the optical observations by John- 
stone, Naylor, & Fabian (1991), who reported that the /-band surface brightness of A1795 
consists of two components (though the extended one is quite elongated), crossing over at 
-Ropt 100/150^ kpc; this is close to i?opt- Furthermore, Mexccss estimated above is similar to 
the virial mass of the cD galaxy (Mviriai — 1.1 x 10^^ ^o), derived from the stellar velocity 
dispersion (~ 300 km sec~^; Oegerle & Hoessel 1991; McNamara & O'Connell 1993). In 
comparison, we have Rx ~ -Ropt ~ 60 h^Q kpc and Mexcess ~ -^viriai ~ lO"*^^ Mq in the Fornax 
cluster. 

The hierarchy interpretation relies much on the fact that the logarithmic curvature of 
Mtot(-R) becomes apparently concave over the region where the shouldcr-likc structure is 
present {R — 150 — 350 kpc). However, this may be a model dependent artifact, produced 
e.g. by cross-over between the two components in the employed double-/? model. As a 
matter of fact, it is rather difficult to assess the reality of this feature, due to complex error 
propagation from the data to the mass curve. Instead of attempting to do so, we here resort 
to fit the observed 0.7-10 keV GIS brightness profile by a rather arbitrary empirical function 
of the form 



where >S'(0), ri, r2, Ci, and C2 are parameters. The best-fit of this empirical model is found 
with n = 7'.72, r2 = 0'.60, ci = 1.67, and C2 = 1.05, yielding xV^^ = 23/15, which falls 
between those from the single- /3 fit and the double-/? fit. The mass profile derived from this 
ad- hoc model, via back projection and equation (4), is also plotted in Figure 8. It reveals 
a very similar feature with concave curvature, even though equation (6) does not involve 
particular angular scales corresponding to 2' — 3' and the logarithmic derivative of equation 
(6) changes monotonically with r. This adds to the reliability of the concave feature seen in 
Figure 8, and reinforce the interpretation of our results in terms of the hierarchical potential 
structure between the cluster and the cD galaxy. 

Finally, some mention may be made on cooling fiows. The dominance of hot ICM in the 
cluster center region, with a smaller amount of cooler component coexisting with it, is in 
itself consistent with the picture of inhomogeneous cooling fiows (Fabian 1994). The short 
coohng time of 2.5 Gyr within the central 80 kpc, as estimated by Edge, Stewart, & Fabian 
(1992), remains essentially unchanged by our results. However, causal relation between the 




(6) 



nested potential structure revealed here and the cooling flows in A1795 remains unclear. 
The hierarchical mass distribution may have formed first during the cluster formation, and 
provided a favorable environment for cooling flows to take place. If so, the cooling flow 
may be regarded as the result, rather than the cause, of the central excess X-ray emission. 
Alternatively, a certain fraction of the excess mass within Rx may have deposited from the 
cooling flow itself. 

In summary, we have discovered that the central excess X-ray emission from A 1795 is seen 
in the hard X-ray band above 3 keV, with a similar fraction as seen in the soft X-rays. This 
property is mainly due to a central excess in the gravitating mass, by Mexcess ~ 3 x lO-*^^ Mq, 
which is hkely to be associated wit the cD galaxy. 
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R. F. Mushotzky, S. Sasaki, and Y. Suto for helpful comments. This research has made use 
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Figure captions 



Figure 1: Background-subtracted GIS2 + GIS3 images of A1795 in (a) 0.7-3 keV and (b) 
3-10 keV. The images have been smoothed with two-dimensional Gaussian filters {a — 0'.5). 
No correction has been made for the XRT vignetting or partial shadows due to the detector 
support ribs. The contours are with the same logarithmic steps in the two images. The SIS 
field of view is indicated by the dashed line in each plot, and the position of the Seyfert 1 
galaxy EXO 1346.2-1-2645 is marked by a cross. 

Figure 2: (a) Spatially sorted GIS2 + GIS3 spectra of A1795 for 0' - 2' and 2' - 12' regions. 
The background has been subtracted, but the detector response has not been removed, (b) 
Same as (a) but for the SISO -|- SISl spectra. 

Figure 3: (a) The GIS spectrum in the r < 2' region (same as shown in Fig. 2a), displayed 
in the form of channel-by-channel ratios to the simulated GIS spectrum. The Montc-Carlo 
simulation was done taking fully into account the PSF effects, assuming a constant ICM 
temperature of 6.2 keV and a constant abundance of 0.3 solar. The model birghtness profile 
was represented by a single /3-model with a core radius of 0.'8 and /3 = 0.56. Eerror bars 
represent ±1 a due to photon counting statistics in the actual spectrum, while those associ- 
ated with the simulated spectrum is negligbly small, (b) Same as (a), but for the r = 6' — 8' 
region. 

Figure 4: (a) Background-subtracted radial count-rate profiles of A1795 obtained with GIS2 
-I- GIS3 in 0.7-3 keV and 3-10 keV. The profiles are corrected for neither the XRT vignetting 
nor the XRT PSF effects. The 0.7-10 keV XRT PSF at the detector center is also plotted 
for comparison, (b) Same as (a) but for SISO. (c) The radial profiles of the 3-10 keV to 
0.7-3 keV brightness ratio (or the spectral hardness ratio) of A1795, derived with G1S2 -|- 
GIS3 and SISO. (d) The same GIS ratio as in (c), but with a 3-10 keV radial profile of a 
point source, whose integrated count is 3% of that of total cluster, added to or subtracted 
from the 3-10 keV GIS data. 

Figure 5: Single-/? fittings to the radial count-rate profiles of A1795 over the whole radius 
range of r < 12'. (a) The ASCA GIS profile in 0.7-3 keV, where crosses represent the data 



while circles show the best-fit model convolved with the XRT + GIS PSF. The lower panel 
shows the fit residuals, (b) The same as panel (a) but, in the 3-10 keV range, (c) The 
ROSATPSPC profile in 0.2-2 keV fitted in the same manner. 



Figure 6: The same as Figure 5, but the fit range is restricted to the outer region of r^5'. 

Figure 7: Two-/3 fittings to the brightness profiles of A1795 over the whole radius range of 
r < 12'. Meanings of panels (a) through (c) are the same as those in Figure 5. Squares 
and triangles show the wider and narrower /3-components, respectively, while circles indicate 
their sum. Panel (d) is the fit to the ROSAT PSPC data, using the X-ray surface brightness 
profile which is predicted by the universal halo model of equation (2) in the text. 

Figure 8: Radial distributions of the integrated total gravitating mass and the ICM mass. 
Cases 1 and 2 are based on the best-fit double-/? model in 0.7-10 keV (Table 4), assuming 
either a constant ICM temperature of 6.2 keV (case 1) or a temperature gradient given by 
equation (5) in the text (case 2). Case 3 is derived by replacing the double- /3 model with 
another analytical surface-brightness model described with equation (6) in the text, with the 
temperature assumed to be constant at 6.2 keV. 



Table 1: Single-temperature (IT) and two-temperature (2T) Raymond-Smith fits to the SIS 
and GIS spectra for the innermost region of 3' radius.* 



Model 




kThot 




Abundance 






(1020cm-2) 


(keV) 


(keV) 


(solar) 




IT 


3t 


6.2 fixed 




0.3t 


645/278 


IT 


3±1 


5.16 ±0.11 




0.34 d= 0.02 


453/277 


2T i 


5±1 


6.5412:1^ 




0.32 ± 0.03 


377/275 



* The SIS and GIS spectra are fitted simultaneously. 

t Errors were not calculated because of too poor a fit. 

t Two-temperature (2T) fits, where the two components are assumed to have the same abundances 
and a common absorption. Abundance rataios are assumed to be solar. Errors represent the single- 
parameter 90% confidence limits. 



Table 2: Single-/? fits to the count-rate profiles in the entire r^l2' range.* 

Detector Energy (keV) Core radius P X^/^ 

ASCA GIS 0.7 - 10 0'.9 0.62 45.5/17 

0.7-3 0'.8 0.60 45.9/17 

3-10 I'.O 0.61 29.6/17 

ROSAT PSPC 0.2-2 O'.S 0.56 49.4/16 



* Confidence ranges of the parameters are not shown, since none of the fits are acceptable. 



Table 3: Single-/? fits to the radial count-rate profiles in the outer region (r^5'), with con- 
servative errors. 



Detector 


Energy (keV) 


core radius 


/3 




central excess* 


ASCA GIS 


0.7-10 


3'.1 ±0'.l 


0.81 ±0.05 


6.4/8 


0.22 




0.7-3 


3'.2±0'.1 


0.79 ±0.05 


6.2/8 


0.28 




3-10 


3'.1 ±0'.l 


0.80 ± 0.05 


2.5/8 


0.20 


ROSAT PSPC 


0.2-2 


3'.8 ± 0'.3 


0.90 ± 0.03 


4.3/9 


0.32 


Einstein IPC^ 


0.5-3 


2.79 ± 0.'93 


0.72 ± 0.08 




0.36 



* Data excess in r^5' above the best-fit model, relative to the total data counts in r < 12'. 
^ Prom Jones & Forman 1984. 



Table 4: Double-/? fits to the r < 12' brightness profile, with conservative errors. 







Compact component 


Extended component 




Detector 


Energy (keV) 


Core radius 




Core radius 






ASCA GIS 


0.7-10 


0'.9±0'.1 


0.83 ±0.05 


3'.2±0'.1 


0.81 ±0.05 


10/14 




0.7-3 


l'.0±OM 


0.90 ±0.05 


3'.3±0'.1 


0.80 ±0.05 


10/14 




3-10 


0'.8 ± OM 


0.79 ± 0.05 


3'.4±0'.1 


0.81 ±0.05 


8/14 


ROSAT PSPC 


0.2-2 


I'.O ± OM 


0.83 ±0.01 


3'.4±0'.1 


0.80 ±0.01 


7/13 



Table 5: Luminosities of the central excess, the cool component, and the total cluster 
emission.* 



Energy (keV) L\^,^^ LI^^^^^ Lexcess/i^totai L 



cool 



0.5-10 17.80 ±0.12 8.37 ±0.54 0.47 ±0.03 1.43 ±0.39 

0.5-3 10.00 ±0.07 4.92 ±0.30 0.49 ±0.03 1.20 ±0.33 

3 - 10 7.80 ±0.05 3.45 ±0.31 0.44 ±0.04 0.23 ± 0.06 

* Luminosities arc derived through the combined fitting of the ASCA GIS and SIS data, in unit of 
^''g s ^ with 90% confidence errors. 
Integrated within 12' of the centroid. 

t Lexcess refers to the luminosity of the narrower component in Table 4. 
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